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Abstract. We investigate theoretically the electromechanical properties of freely 
suspended nanowires that are in tunnelling contact with the tip of a scanning tunnelling 
microscope (STM) and two supporting metallic leads. The aim of our analysis is 
' to characterize the fluctuations of the dynamical variables of the nanowire when a 

temperature drop is mantained between the STM tip and the leads, which are all 
■<^J- | assumed to be electrically grounded. By solving a quantum master equation that 

describes the coupled dynamics of electronic and mechanical degrees of freedom we find 
that the stationary state of the mechanical oscillator has a Gaussian character, but 
f — ' that the amplitude of its root-mean square center-of-mass fluctuations is smaller than 

would be expected if the system were coupled only to the leads at thermal equilibrium. 
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1. Introduction 

The possibility to detect and control the motion of nanometer-sized-mechanical 
resonators by coupling them to mesoscopic electronic devices has generated a 
considerable research effort in recent years [TJ. In particular, the possibility to 
use such nanoelectromechanical systems (NEMS) as ultrasensitive sensors of, for 
example, displacement and mass, have been demonstrated in a number of works [2j[3]. 
Independently of the specific type of electronic device considered in the different schemes, 
a common feature that has emerged from these studies is that the electronic subsystem 
must be out of thermodynamic equilibrium in order to function as an ultrasensitive 
measurement tool in combination with the mechanical subsystem. This observation 
naturally raises the question of how the dynamics of the mechanical subsystem is affected 
by the nonequilibrium environment created by the mesoscopic electronic device to which 
it is coupled. 

It is known from statistical mechanics that the displacement and momentum 
fluctuations of a quantum harmonic oscillator (the basic model for any movable structure 
that could be included in a NEMS), which is coupled to a thermal bath in equilibrium 
at temperature T are described by Gaussian distribution functions, whose widths are 
proportional to coth 1//2 [hioj '2k gT) where u is the frequency of the oscillator. From 
this formula it follows that the fluctuations have a thermodynamic origin in the high- 
temperature limit, where they are fully defined by the temperature of the thermal bath. 
On the other hand, at low temperatures &bT <C Hoj the fluctuations are completely 
quantum mechanical in nature. 

What kind of changes from this picture could be expected if the oscillator is 
coupled to a nonequilibrium environment? In spite of the difficulties related to the 
definition of temperature for systems that are out of thermodynamic equilibrium, 
several theoretical works show that nonequilibrium fluctuations in the properties of 
nanomechanical oscillators coupled to mesoscopic electronic systems (such as a single- 
electron transistor or a superconducting Cooper pair box) are, to a good approximation, 
still described by Gaussian distribution functions [4, [Sj. 

A remarkable difference between a passive (thermodynamic) environment and an 
active (nonequilibrium) one is that in the latter case the amplitude of the fluctuations 
can be controlled through parameters that characterize the state of the electronic 
subsystem. This external control introduced by the coupling between mechanical and 
electronic degrees of freedom opens the way for the possibility to reach the quantum 
limit of fluctuations even when the temperature is high on the scale defined by the 
quantum of mechanical energy, i.e. hou/hs- 

In the last few years considerable efforts have been made in order to, develop efficient 
procedures to effectively "cool" down the motion of nanomechanical resonators below 
the threshold defined by thermal fluctuations. Most of the proposed schemes strive to 
reproduce the effects of laser cooling of atoms and molecules by purely electronic means 
[3 El [TJ [9] . The general strategy underlying these approaches is based on the coherent 
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control of resonant, energy-conserving transitions between discrete electronic levels. 

Recently we have proposed a fundamentally new scheme to cooling the vibrations 
of a suspended-nanowire based mechanical oscillator (TO]- In contrast to the 
aforementioned cooling schemes, our proposal has the advantage that it does not rely 
on the energy conservation constraint. In particular, we considered a suspended carbon 
nanotube in tunnelling contact with the voltage-biased tip of a scanning tunnelling 
microscope (STM) and the metallic leads at which its ends are clamped. Our analysis 
showed that the average number of quantized vibrational excitations, i.e. vibrons, 
(which is proportional to the root-mean square fluctuations of the center-of-mass 
position) can be reduced by varying the bias voltage within a range of values for which 
the probability for absorbing vibrons during inelastic electron tunnelling processes is 
significantly enhanced over the probability for vibron emission. 

In this paper we demonstrate that the cooling mechanism suggested in Ref. 10 
can work by exploiting the temperature gradient, rather than the potential drop across 
the system. We show that the amplitude of the root-mean square fluctuations of the 
nanotube center-of-mass position is smaller than what would be in presence of only an 
equilibrium thermal bath. This partial suppression of the stochastic fluctuations of the 
nanotube displacements can be interpreted as an effective cooling of the mechanical 
degrees of freedom of the system. Moreover, we have found that this effective cooling 
phenomenon involves simultaneously several low-frequency vibrational modes and not 
only the fundamental one. 

2. Model Hamiltonian 

To be specific, we consider the system sketched in Fig. HJ where a carbon nanotube is 
suspended over a trench between two metallic leads. Its ends are both clamped, while 
a third electrode in the form of the tip a scanning tunnelling microscope (STM) is 
positioned above the nanotube. The suspended carbon nanotube can be considered as 
a quantum dot that is coupled to the surrounding electrodes through tunnel junctions. 
Low-temperature tunnelling spectroscopy studies on freely hanging carbon nanotubes 
have shown that inelastic electron tunnelling can create a non-thermal equilibrium 
population of the vibronic states of the nanotube [TT] . 

In order to analyze the dynamics of the nanotube deflections and the behavior of 
the electronic subsystem in the quantum regime we introduce a model Hamiltonian, 
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Figure 1. Sketch of the model system considered. A carbon nanotube is suspended 
over a trench between metallic leads, while an STM tip is placed a distance d above 
the nanotube. Both the STM and the leads are grounded, while their temperatures 
T s and T L arc different and held constant (T s > T L is assumed). 



H = H e + H m + Ht + He, where the different contributions are given by: 
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a g \a and are annihilation (creation) operators for 



In Eqs. (fTij), (flTj) and (fid 

electrons in the STM tip (a = S), in the leads (a = L) and in the nanotube, respectively. 
The term H e in Eq. ( JTa|) describes the electronic states in the STM tip, the leads (which 
are treated as reservoirs of non-interacting quasiparticles) and in the nanotube. We 
assume that, in the range of temperatures that are relevant for our considerations, only 
one quantized electronic level in the nanotube is involved in the exchange of charge 
with the reservoirs and we denote its energy as E . Such a condition is satisfied if the 
temperatures of the reservoirs are significantly lower than the characteristic difference in 
energy between the quantized electronic levels of the nanotube, which can be estimated 
as A ~ hvp/2£ w 1.7 meV//xm, where Vp is Fermi velocity and i is the length of the 
nanotube. 

The term H m in Eq. ( Ilbl) describes the mechanical degrees of freedom of the 
nanotube. The quantum field u(z) gives the nanotube deflection from the straight 
configuration at point z (that is, the nanotube axis, see Fig. [Q, while tt(z) is the 
momentum linear density and the symbol ' denotes derivation with respect to the 
coordinate z. The displacement and momentum density fields are canonically conjugated 
dynamical variables, that is they obey the commutation relation: [u(zi),fr(zj)} = 
ih5(zi — Zj). 

The parameter p represents the linear mass density of the nanotube, k its bending 
rigidity, and L is the length of the suspended part (notice that L ^ £, where the latter 
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is the total length of the nanotube). The clamping of both the ends of the nanotube to 
the leads can be expressed through the boundary conditions u(±L/2) = u'(±L/2) = 0. 

The tunnelling of electrons through the STM tip-nanotube and the nanotube-leads 
junctions is described by the Hamiltonian operator H T , presented in Eq. ( flc|) . We denote 
by zq the point along the nanotube axis above which the STM is positioned. Both the 
tunnelling amplitudes ts[u(z )] and are assumed to be independent of the electronic 
energy, whereas only t s [u(z )] is a function of the nanotube deflection, as a consequence 
of its dependence on the overlap of the electronic wavefunctions in the STM tip and 
the nanotube. We model this deflection dependence of the probability amplitude of 
tunnelling between the STM tip and the nanotube as ts[u(zo)] = ts exp[u(z )/ X], where 
A is the characteristic tunnelling length of the junction (A ~ 10~ 10 m). 

The effect that the nanotube displacement has on the width of the tunnel barrier 
provides a mechanism of coupling the electronic and the mechanical degrees of freedom 
of the system. In the following, we will refer to that as tunnelling electromechanical 
coupling. 

The last term in the Hamiltonian, Hq, shown in Eq. flldj) . describes the electrostatic 
interaction in the system. Since both the STM tip and the substrate are grounded, an 
electron occupying the state inside nanotube induces a polarization charge of opposite 
sign in the STM tip and hence generates an electrostatic force 3 acting on the nanotube. 
This force can be thought as applied to the point zq and directed towards STM tip (this 
"strongly localized" form of the electrostatic force is a valid approximation if the effective 
radius of the STM tip is negligible with respect to the length of the nanotube). 

The electrostatic interaction described by the operator H c provides another 
mechanism that couples the dynamics of the mechanical and electronic degrees of 
freedom, which we will refer to hereafter as the polaronic electromechanical coupling, 
because of the formal analogy with the interaction term in Hamiltonian of the polaron 
problem. Under the assumptions of uniform charge distribution inside the charged 
nanotube, if the equilibrium distance between the nanotube and the STM tip d is much 
less than the effective radius of the tip, R, then ^ can be approximated by the following 
expression: 



where (3 is a numerical factor of the order of one and e$ the vacuum permittivity. 
From the considerations presented above, it follows that the electronic and mechanical 
subsystems interact through two independent coupling mechanisms, that is the 
tunnelling and the polaronic. The former one results in the change of the nanotube 
momentum on the value h/X when one electron tunnels from the STM to the nanotube 
or in the opposite direction. The second one accounts the difference between the 
equilibrium configurations of the charged and neutral nanotube. Working incoherently 
these two mechanisms of electromechanical coupling give rise to stochastic fluctuations 
of the nanotube center-of-mass position. 

However, as it was shown in Ref. 10, the quantum interplay between them may 
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significantly modify their mutual performance. In order to analyze the consequences 
of the combination of the two coupling mechanisms, it is convenient to introduce the 
eigenmode representation for the nanotube displacement and momentum density fields, 
that is defined by the operators 

1 f L/2 

X n = —j= \ dz<p n (z)u(z) (3a) 

V L J -L/2 

r~ f L/2 

P n = VL dzip n (z)7i(z). (3b) 

J -L/2 

The operators presented in Eqs. f l3a|) and (I3bj) satisfy canonical commutation relations, 
that is [X n , Pi] = ih5 n) i and the complete set of orthonormal functions (p n (z) is given 
by the eigenfunctions of the operator d 4 /dz 4 with the boundary conditions ip n (±L/2) = 
(p' n (±L/2). In this representation the terms H m , Ht and He in the Hamiltonian assumes 
the form: 

H m + H C = J2 (l^Pn + ^Xn) ~ ^cJ2?n(zo)X n , (4) 
n ^ * n 

Ht = J2 [tse E ^ Zo) * n/X Ja qiS + t L al >L c] + H.c, (5) 
<? 

where M = pL is the mass of the suspended part and u n is the eigenfrequency 
of the n-th bending mode. The polaronic term can be removed from the 
Hamiltonian by a suitable unitary transformation, H — > H = UHW, where 
U = exp[ih~ 1 Qc^cJ2 n i Pri(zQ)P n /2Mu 2 ]. However, as additional consequence of this 
transformation, the tunnelling amplitudes turn out to be dependent on the momentum 
operators. The transformed tunnelling Hamiltonian reads 



E 



+ H.c, (6) 



The physical analysis of Eq. (|6]) becomes more transparent after introducing the creation 
(annihilation) operators b' n (b n ) for the elementary mechanical excitations (vibrons) of 
the n-th ben ding mod e: X n = (6+ + b n )x n /V2, P n = ^(v^Xn) -1 ^ - K), where 
Xn/y/2 = \/fr/2Muj n is the position uncertainty in the vibrational ground state of 
the oscillator associated to the n-th mode. In this representation, the part of the 
Hamiltonian which describes the electron tunnelling processes between the STM tip 
and nanotube assume the form: 

t s {e£^ibl+A; bn ) c t aqiS + e E n (^bi+AiK) a tA ; (7) 

q 

where the parameters A^, which characterize the rates of the inelastic electronic 
transitions with absorption and emission of the vibronic quanta, are given by: 
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From the expressions of the parameters shown in Eq. (JSJ), it follows that the vibron 
emission processes are suppressed with respect to the absorption ones when electrons 
tunnel from the STM tip to the nanotube, whereas the absorption is suppressed and the 
emission promoted during the transitions in the opposite direction (which are described 
by the rightmost term of the transformed tunnelling Hamiltonian shown in Eq. (JTj)). 

Moreover, one can achieve complete suppression of the vibron emission (absorption) 
for the given mode by varying the value of the electrostatic force which, according 
to Eq. (j2J), is controlled by the equilibrium distance between the STM tip and the 
nanotube, d. If after tunnelling from the STM tip, the electrons tunnel immediately 
off to the leads, the rate at which vibron emission processes occur will be substantially 
reduced and therefore, in the stationary regime, one can expect that the number of 
vibrons will be close to zero. 

However, in order to drive a certain vibrational mode to its quantum ground state 
(which corresponds to an average number of vibrons much smaller than 1), its frequency 
must satisfy the condition ui n = ^/2X^s/h. In the stationary regime, the root-mean- 
squared deviation of the center-of-mass position (which expresses the fluctuations of the 
mechanical state of the nanotube around its equilibrium configuration) is given by the 
square root of the following expression: 



Taking into account that for the doubly clamped nanotube wo/w„ ~ (n + l) -2 and that 
V 9 2n+i(0) = Eq. dHJ) indicates that the of center-of-mass fluctuations are mainly defined 
by the average number of the vibrons in the fundamental mode. Therefore, in order to 
suppress such fluctuations, this number shoud have the minimum possible value. On 
the basis of these considerations, in the rest of this paper we restrict our attention to the 
fluctuations of the fundamental bending mode, which will be considered as a quantum 
harmonic oscillator. 

3. Quantum master equation 

In order to perform a quantitative analysis of the phenomena discussed above, we 
start from the Lioville-von Neumann equation for the density matrix operator, which 
represents the state of the whole system 



In Born approximation with respect to the tunnelling amplitudes tg, ti, Eq. (JTDjl can 
be recast in the following integral form: 




(9) 




(10) 



dp 1 



[#o,p]-4 / dt'[H T (t),[H T (t'),p(t')}} 



(11) 



dt ih 



J — oo 
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where H = H e + H m + H c , and A(t) = e iH t/h^ e -iH t/h_ Taking into account that 
the coupling between the nanotube and the electronic reservoirs (that is, the STM 
tip and the leads) is weak enough so that any back-action of the nanotube on their 
physical states is negligible, we can use the Ansatz: pit) ~ a(t) ®Ps®Pl- The operator 
a(t) is the reduced density matrix operator, which is defined as a(t) = Trg+L[p(t)] and 
represents the electronic and mechanical state of the oscillator. The density matrices p a 
describe the STM tip (a = S) and the leads (a = L) as electronic reservoirs at thermal 
equilibrium with temperatures T$ and Tl, which means that: 

^{al, a a q , a pa) = (l + exp{{E qta) - p a )/k B T a }) 1 = f a (E qt0l - p) 

where the p is the chemical potential of the STM tip and the leads, which is supposed 
to be the same. 

In order to describe the nanotube dynamics it is convenient to project the reduced 
density matrix onto the subspaces corresponding to the electronic level in the nanotube 
being occupied or unoccupied. That amounts to multiply the operators c^c and cc^ 
to Eq. (TTTil and trace over the electronic degrees of freedom of the nanotube. After 
this procedure, we obtain two coupled equations for the operators a = Tr c (acc^) and 
o\ = Tr e (o"c^c) which, in the high-temperature limit hw/ksTL <C 1 turn out to be 
local in time |12j . For small displacements of the nanotube around the equilibrium 
configuration, the tunnelling amplitude ts[X) can be linearized, so that the equations 
of motion for cr ,i have the form: 

^ = -~[H m ,vi] + §[X,°i\ ~ I>i " T 5 U + j{a h X} + ^{<xi,A 2 }^ + 
+ T+ (v + j{a , X} + t^Wo, A 2 } + -^Acr A^ + T+^o + Cy<n 

d ~2 = ~\[HmM + rz^i + r s L + \{a u X} + ^Wi,X 2 } + ^ 2 Xa x X 

- T% (c + i{a , A} + ^{<x , X 2 }^j - T+ao + £ 7 <x , (12b) 

where T+ = T a f a (E ), T~ = T a [l — f a (E )} and T a = 2Tih~ 1 \t a \ 2 v a , v a being the density 
of states at the Fermi energy in the STM tip (a = S) and in the leads (a = L). 

The operator £ 7 in Eqs. (11 2 aft and ( 112bl) models the relaxation of the oscillator 
towards thermal equilibrium with the phononic bath in the leads, a process characterized 
by the rate 7 = ojq/Q, where Q is the quality factor. On the basis of general 
considerations regarding quantum dissipative systems p2], £ 7 can be explicitly written 
as: 

£>] = -?L[X,{P,o-}] -^Lcoth(fko /2k B T L )[X, [A, a]]. (13) 

The mechanical state of the suspended nanotube can be characterized through the 
operator cr + = a + ox, whose evolution in time is fully determined by Eq. (112aj) 
and (112bl) . once that the operator cr_ = a — (j\ is introduced. Furthermore, in order 



(12a) 
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to describe the stationary state of the oscillator, it is convenient to introduce the 
dimensionless operators x = X/xo, P = XoP/h an d express the operators a± in the 
"Wigner function representation" [H], which is defined as: 

W ± (x,p,t) = / -i e -«W(x-$|a ± (t)|a; + (14) 

From Eqs. (Il2ap and (112bl) . it follows that the Wigner functions corresponding to the 
stationary solutions of the quantum master equations for the operators a + , cr_ satisfy 
the equations: 



(p8, - xd r )W + = |^3 P W + + |a,TF_ + e£±($W + + ^%W-+ 



+ ^(PW + ) + + <KQ-',4)W- (15a) 

(pd x - xd p )W. = e p —^d p W. + -£d p W + - (Ar s + 2e t AT s x) W+- 

ZL s Z 

- (r s + 2e t T s x) + 0(Q-\ e\)W. (15b) 

where e t = (p (z )xo/ \, £ p = (Po(z )Q/xoMcOq, while the parameters related to the 
tunnelling processes have been rescaled in units of coq, that is AT a = (r+ — T~)/uo, 
where a = S,L, Ar s = (Ar 5 + AT L )/u and T s = (r s + T L )/u . 

Eqs. (I15aj) and (115bl) can be solved by means of a perturbative expansion in the 
small coupling constants e p ~ e ( < 1 and inverse quality factor Q^ 1 ~ e^, e|. At the 
zero-th order in the small parameters, the solution of Eqs. (I15a[) and (|15ap . has the form 
W, = w(A), W _ = — (Ar 2 /r s )w(y4), where A = a/x 2 + and w is an arbitrary 
function. 

The necessary and sufficient condition for w(A) to be a good zero-order 
approximation of W + is that any deviation from w(A) is at most of order (ef, e 2 , Q~ l ) in 
comparison to w(A). Then, by replacing the expressions of W± up to the second-order 
corrections in the small parameter into Eqs. fl 1 5 aj) and (115bj) . and neglecting all the 
contributions except the zero-th order ones, we obtain an equation that determines the 
form of w(A). 

The first order corrections to the function = w(A) vanish after performing the 
transformation x — > x — x (where x = (e p /2)(l — AT^/T^)), which can be considered 
as a shift of the reference frame that accounts for the nanotube deformation induced by 
the average electrostatic force. 

4. Results 

Following the perturbative procedure described above, it turns out that, for the function 
w(A) to be an appropriate approximation for W + , it must satisfy the following first- 
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order linear differential equation: 

*;r E (i - (Ar s /r s) 2 ) + 4r s r s -Ar s Ar s + i / 



1 + 4 T s 2Q \2k B T L 



e p e t r L ATs-T s Ar L , 1 



- 2 (It tct + q) Am < 16 > 

Now we focus our attention on the regime in which the transport of charge from the 
STM tip to the leads is activated by the gradient of temperature between the tip and 
the leads. In order to have an appreciable rate of tunnelling between the STM tip and 
the nanotube, the STM tip temperature must satisfy k B T s ~ E Q , so that the electronic 
states at the energy of the electronic level of the nanotube have a good chance to be 
populated, that is fs(E ) ~ 1/2. At the same time, the temperature in the leads 
should be much lower than T s , so that the electronic states in the leads at the energy 
corresponding to the electronic level in the nanotube have a good chance to be empty, 
which means Jl(Eq) ~ 0. Furthermore, the temperature of the leads should be high on 
the scale defined by the vibrational quantum energy, that is k B T^ 3> frwo. 

For what concerns the best cooling performance, we already observed that it is 
expected to be achieved if the electrons, after having tunneled to the nanotube from the 
STM tip, tunnel quickly to the leads rather than being transferred back to the STM 
tip. That corresponds to the situation in which Ts/Tl <C 1. Under these conditions, 
the quasi-distribution function w(A) has the following Gaussian form: 

w(A) = —e-£, (17a) 



2(i+f|) 1 r s r L 



+ 



+ + C0 \h{nw Q /2k B T L )^- L ^ 



(17b) 



2(1+1*) 



We remark that, by virtue of the Gaussian form of the stationary state, the root-mean- 
square fluctuations of the nanotube center-of-mass position are just proportional to the 
width of the quasi-distribution function, i. e. (X 2 ) 1 / 2 = 9/y/2. 

For an oscillator coupled to an equilibrium environment, in the high-temperature 
limit, 6 2 becomes proportional to the temperature of the phononic bath in the leads, 
9 2 ~ Tl (in agreement with Einstein's relation), while it reduces to 1 for the w(A) that 
describes the quantum fluctuations of the oscillator in the ground state. 

The denominator of the ratio that defines 8~ 2 describes the diffusion in the energy 
space of the oscillator induced by the stochastic electronic tunnelling processes, whereas 
the numenator represents the effective damping generated by them (see Eq. (|17bl) ). 

In order to understand the physical origin of the stationary state of the oscillator 
determined by the temperature-activated electron tunnelling, let us consider the limit 
Q — > oo, that is the situation in which the nanotube is decoupled from the equilibrium 
environment. In this limit the quasi-distribution function is determined only by the 
rates of the inelastic tunnelling transitions induced by the polaronic and tunnelling 
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electromechanical couplings, which are characterized by the parameters e p and e t , 
respectively. The stationary state in this limit is characterized by a width given by 



where 77 = ^J2e 2 /e 2 {l + It follows from Eq. ( fTHI) that for both strong and weak 
polaronic coupling 9 is larger than 1, which corresponds to a state that is far from 
the quantum ground state. Nevertheless, in a suitable range of values of e p , 9 2 is 
significantly smaller than coth(huo / 2k bT^), which means that the stationary state can 
be interpreted as a thermal state characterized by an effective temperature smaller than 
T L . Correspondingly, the root-mean-square fluctuations of the nanotube center-of-mass 
position are smaller than the value determined in the thermal equilibrium case, therefore 
the mechanical subsystem is effectively cooled. The minimum value that 9 2 can reach as a 
result of the interaction of the nanotube with the nonequilibrium electronic environment 
is given by 6 2 min = ^2(1 + T|) > 1. 

From Eq. ( IT8l) . one can see that the effective cooling of the mechanical 
vibrations induced by the nonequilibrium environment requires the presence of both the 
mechanisms of electromechanical coupling in order to exist. In the case in which only 
one mechanism is active, it follows from Eq. (117b|) that it can only generate diffusion in 



energy space, which results in a broadening of the quasi-distribution function. We stress 
that this behavior is characteristic of the nonequilibrium situation considered here, since 
in the case of coupling with an equilibrium environment, both damping and diffusion 
are present, independently of the type of interaction. 

The stationary quasi-distribution function generated by the temperature drop from 
the STM tip to the leads is plotted with respect to the amplitude A (expressed in units of 
Xo) i n Fig- El for different values of the quality factor. Furthermore, we can compare the 
size of the root-mean-square fluctuations of the nanotube center-of-mass position in the 
ground state, the thermal equilibrium state and the electronically-induced stationary 
state as a function of the equilibrium distance between the nanotube and the STM tip, 
d, for different quality factors, as shown in Fig. [3l The curves shown in Figs. (|2]) and ([3]) 
indicates that the interaction between the electron tunnelling current and the oscillator 
can be interpreted as an effective cooling of the mechanical degrees of freedom. 

In order to detect experimentally the cooling effect predicted above, the most 
direct approach consists in the measurement of the root-mean-square fluctuations of 
the nanotube center-of-mass position. Regarding this point, it is has been argued since 
a long time ago that the STM (combined with a current amplifier) can provide the 
basic building block for a quantum-limited position displacement sensor. The tunnelling 
current that can be measured at the output of such a device contains information about 
the displacement of the mechanical system under investigation but, at the same time, 
perturbs it with a very small back-action force, being this mainly due to the random 
momentum transfer associated with the tunnelling electrons [15J. 
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Figure 2. Comparison between the Wigner functions corresponding to the ground 
state of the oscillator, the thermal equilibrium state at temperature Tl and the 
stationary state induced by the nonequilibrium electronic environment for different 
quality factors. Values of the relevant parameters: u>q = 10 9 Hz, T$ — 5 • 10 6 Hz, 
T L = 10 8 Hz, e t = 0.19, e p = 0.14. 




Figure 3. Comparison between the fluctuations of the root-mean-square fluctuations 
of the center-of-mass nanotube position y/ (X 2 ) calculated with the Wigner functions 
corresponding to the ground state of the oscillator, the thermal equilibrium state at 
temperature Tl and the stationary state induced by the nonequilibrium electronic 
environment as a function of the equilibrium distance between the STM tip and the 
nanotube, d, for different quality factors. Values of the relevant parameters: ui = 10 9 
Hz, T s = 10 8 Hz, T L = 10 8 Hz, A = 10- 10 m, e t = 0.19. 
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In conclusion, we have studied the coupled dynamics of the mechanical and 
electronic degrees of freedom of a suspended-nanowire-based NEMS wherein the movable 
element is in tunnelling contact with the tip of an STM and two supporting metallic 
leads. Our analysis shows that, in the regime in which the electron transport is 
activated only by the temperature difference between different parts of the device, an 
effective cooling of the mechanical degrees of freedom can be achieved. This result 
depends crucially on the interplay between tunnelling and electrostatic coupling that 
characterizes the system considered here. The interaction of the suspended nanowire 
with the nonequilibrium environment provided by the tunnelling current reduces the 
amplitude of the root-mean-square fluctuations of the center-of-mass position of the 
nanowire. This effect could be in principle detected experimentally thanks to the almost- 
quantum limited sensitivity of the STM as a displacement sensor. 
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